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ABSTRACT 

The  present  article  is  to  examine  the  joint  injluence  of  internal  heating  and  variable  viscosity  on  the  onset  of  double 
diffusive  convective  motion  using  perturbation  technique.  Linear  stability  analysis  is  performed  and  it  is  assumed  that 
the  upper  surface  of  a  fluid  layer  is  deformably  free  and  that  viscosity  dependency  is  believed  to  be  exponential.  The 
boundaries  are  known  to  be  rigid,  but  permeable,  and  insulated  to  jluctuations  in  temperature.  From  results  of  an 
increasing  the  viscosity  parameter,  the  system  shows  destabilizing  effect  and  the  system  will  be  stabiliied  by  that  internal 
heat  source  power.  It  is  also  revealed  the  injluence  of  double  diffusive  coefjicients.  The  effect  of  the  thermal  diffusion  is 
found  to  have  a  destabilizing  reaction  on  the  system,  whereas  the  opposite  reaction  is  noted  with  an  increase  of  thermo- 
dijjusion  parameter. 
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1.  INTRODUCTION 

Dufour  diffusion,  also  known  as  thermo-diffusion,  and  Soret  diffusion,  also  known  as  thermal  diffusion,  is 
significant  in  both  non-Newtonian  and  Newtonian  heat  convection  and  mass  shifts  and  are  often  found  in  high- 
speed  aerodynamics  and  chemical  process  engineering.  Such  findings  are  important  in  combining  heat  and  mass 
transfer  for  the  intermediate  molecular  weight  gasses  in  binary's.  Hurle  and  Jakeman  (1971)  demonstrated  the 
thermosolutal  convection  induced  by  Soret  using  a  methanol-water  mixture,  both  experimentally  and  theoretically. 
Instead,  using  a  salt  solution,  Caldwell  (1970)  extended  the  study,  while  Platten  and  Chavepeyer  (1973)  continued 
previous  studies  using  a  combination  of  water-ethanol.  The  linear  stability  of  experimental  thermal  diffusion 
convection  in  an  ethanol-water  mixture  was  investigated  by  Knobloch  and  Moore  (1988)  under  various  boundary 
conditions  with  a  focus  on  the  Biot  number.  Taking  into  account  the  Soret  interaction  with  other  physical  causes, 
the  thermo  capillary  instability  in  a  binary  tluid  at  the  start  of  convection  was  studied.  (Bergeon  et  al.  (1988), 
Slavtchev  et  al.(  1 99 1 ).  Saravanan  and  Sivakumar  (2009)).  Chand  and  Rana  (2015)  recently  studied  the  Soret  effect 
in  the  presence  of  a  vertical  magnetic  field  in  a  nanotluid  layer. 

Some  researchers  assume  that  fluid  has  a  constant  viscosity  or  may  possess  temperature-dependent 
viscosity  (viscosity  decreases  exponentially  with  temperature)  which  may  affect  the  stability  convection.  Palm 
(1960)  initiated  the  study  of  variable  viscosity  on  a  steady  convection.  Other  researchers  also  studied  the  variable 
viscosity  effects  in  different  problems  where  Torrance  and  Turcotte  (1971)  and  Stengel  et  al.  (1982)  studied  in 
Benard  instabilities  and  Slavtchev  et  al.  (1999),  Cloot  and  Lebon  (1985)  and  Kozhoukharova  and  Roze  (1999)  in 
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Marangoni  instabilities.  Cloot  and  Lebon  (1985)  studied  the  steady  Marangoni  convection  with  undeformable  surface  and 
in  microgravity.  Abidin  et  al.  (2008)  and  Arifin  and  Abidin  (2009a,2009b)  studied  the  temperature-dependent  viscosity 
effect  together  with  others  effects  such  as  the  feedback  control,  deformable  surlace  and  boundary  effect  in  a  fluid  layer. 
White  (1998)  did  an  overall  study  where  he  studied  both  theoretical  and  experimental  in  a  Benard  convection  meanwhile 
Manga  et  al.  (2001)  compared  the  temperature  effect  experimentally  with  boundary  layer  models. 

Temperature  dependent  viscosity  also  has  been  integrated  into  other  convection  system  where  Franchi  and 
Straughan  (1992)  included  the  effect  in  a  micropolar  uid  and  Ramirez  and  Saez  (1990)  in  a  porous  medium.  It  stated  that  a 
higher  temperature  in  both  systems  will  cause  the  critical  Rayleigh  number  to  decrease.  Ramirez  and  Saez  (1990)  stated 
that  temperature-dependent  viscosity  should  be  taken  into  account  for  every  case  studied  since  the  effect  has  a  huge  impact 
on  the  instability  of  convection.  A  similar  result  obtained  by  Lu  and  Chen  (1995)  where  the  stability  was  enhanced  by 
decreasing  the  temperature. 

Friedrich  and  Rudraiah  (1984)  have  shown  that  the  convection  can  be  suppressed  with  a  sufficient  non-uniform 
temperature  gradient.  Bhattacharyya  and  Jena  (1984)  found  a  destabilizing  effect  of  the  heat  source.  Char  and  Chiang 
(1994)  showed  that  internal  heat  generation  has  a  significant  impact  on  the  stability  of  Benard-Marangoni  convection. 
Through  a  rigorous  investigation,  they  found  that  the  system  is  destabilized  by  an  increase  in  internal  heat  generation. 
Capone  et  al.(201 1)  investigated  the  impact  vertical  through  flow  on  the  onset  of  double-diffusive  penetrative  convection. 
In  a  rotating  anisotropic  porous  sheet,  Bhadauria  et  al.  (2011)  studied  the  influence  of  internal  heat  generation  on  the  onset 
of  natural  convection  using  a  weak  nonlinear  model.  Very  recently,  Yadav  et  al.  (2012,2015)  and  Wakif  et  al.{ 2016) 
examined  the  effect  of  internal  heat  source  in  a  nanofluid  layer. 

In  a  fluid  layer  with  internal  heating,  double-diffusive  convection  despite  its  value  in  many  hands-on  industries. 
The  current  findings  attempt  the  effects  of  internally  heating  fluid  layer  in  the  presence  of  thermo  and  thermal  diffusion. 
We  perform  a  stability  analysis,  and  the  obtained  eigen  value  problem  is  solved  using  the  perturbation  procedure. 

2.  CONCEPTUAL  MODEL 

Two  horizontal  layers  of  quiescent  double  diffusive  binary  fluid  with  thickness  d  is  heated  from  below  where  the 
temperature  difference  is  represented  by  A T .  The  upper  free  surface  of  the  fluid  layer  with  its  location  being  free  of 
deformities  z  —  (I  +  Q.  [  X,  V,  t)  and  the  z-axis  pointing  vertically  upwards  opposite  the  direction  of  gravity.  For  a 
Boussinesq  approximation,  we  assumed  the  constant  physical  fluid  properties  except  the  density,  p  and  surface  tension, 

<J  (for  Marangoni  convection)  to  vary  upon  temperature,  T  and  solute  concentration,  S  Chen  and  Su(1992)  and  takes 
the  form 

&  =  <?0-<?,(T~T0)  +  cts(S^S0)  (1) 

P  =  P  o[}  -  P ,  (T  -  T0)+  P  s  (s  -  S  (2) 

Here,  £T0  and  p{]  are  the  values  at  the  reference  concentration,  S{)  and  the  reference  temperature,  T0  .  (7t  and 
pt  are  the  rate  of  change  with  temperature  and  <7S  and  ps  are  the  rate  of  change  of  density  with  concentration.  Due  to  the 
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temperature-dependent  viscosity  in  a  binary  mixture,  the  kinematic  viscosity,  [1  follows  Hilt  et  al.  (2014)  in  the  form 


=  exp[-A(r-r0)]  (3) 

and  fl{)  is  the  dynamic  viscosity  corresponding  to  a  temperature  equal  to  the  mean  of  temperature  at  the 
boundaries. 


Free  boundarv 


z  —  d 


z=d  +  Ofx.y.r) 

Fluid  layer 


z  ♦ 
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Figure  1:  Physical  Configuration. 

3.  MATHEMATICAL  FORMULATION 

The  governing  equations  of  motion  are: 

vv  =  o 


Po 


c)V 

dt 


-  + 


(v-v)v 


~Vp  +  p0g  +  2V. 


//(v.v  +  v.vr) 


dT 

dt 


dS 


+(v -v)t  =  KV2T + DtsV2T + q 


^  +(V-V)s  =  KsV2S  +  DSTV2T  +  q 


(4) 


(5) 


(6) 


(7) 


where  the  variables  are  represented  as  follows;  V  =  (w,V,  w)  is  the  velocity,  S  (solute  concentration),  p 

(pressure),  k  (thermal  diffusivity),  ks  (solutal  diffusivity),  DTS  (Soret  diffusivity)  and  lastly,  Dsr  is  the  Dufour 
diffusivity.  The  basic  state  of  the  fluid  is 

(u,v,w,p,T,p,fi,S)  =  [0,0,0,pb{z),Tb(z),pb{z),ph{z),Sb(z)']  (8) 

Infinitesimal  disturbances  are  implemented  to  test  the  stability  of  the  basic  solution 

(u,v,w,p,T,p,ju,S)  =  [0,0,0,ph(z),Tb(z),pb{z),Jub{z),Sb  (z)\  +  \u  ,v  ,w ,  p  ,T' ,  p,p,  5’]  (9) 


And  using  equations  (8)  and  (9),  we  obtain 
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1  d 


V2w  =  fV4w+2^V2^  +  ^ff(V2w-2V2w)  +  RV:T  +  RsV2hS 
pr  dt  dz  dz  dz  v  7 


^  =  VfT  +  Df  V\S  +  w[  1  -  Ns(l  -  2 Z)\ 
at 


=  SrV2hT  +  LeV2hS  +  w 
ot 


(10) 


(11) 


(12) 


D  ag  ATd'  psg  ASd3 

Where  R  = -  (Rayleigh  number),  Rs  = -  (Solutal  Rayleigh)  number. 


V  K 


VK 

K, 


dst  a  s 


Ns  =  qd  /2k(Tq—Tu)  (dimensionless  heat  source  strength),  Le  =  — —  (Lewis  number),  Sr  = 


(Soret 


number)  and  DF  = — (Dufour  number)  and  Pr  =  —  (Prandtl  number).  V1  =-V2h +  d2  /  dz2  is  the 
Laplacian  operator  with  V),  =  r) 2  /  dx2  +  d2  /  dy2 .  The  function  /  representing  variable  viscosity  function,  is  defined  as 


(  iY 

( v  ] 

B 

z  — 

,  B  = 

V  ^J. 

(.  Vrmn  / 

/  =  exp 


The  respective  boundary  conditions  are 

=  0,  w  =  0  and  =  0  at  z  =  0 
oz  OZ 

dQ  w,^  +  Bi[T-(Ns  +  l)nl  =  0 

I—  —I 


dt  dz 


f 


(  wi  \ 

— _ V2 

9z2  * 


w 


=  MV2[r-(l  +  M5)^] 


at  z  =  1 


at  z  =  1 


(13) 


(14) 


(15) 


(16) 


fCr 


Pr  dt 


(  ^2 


dz 


+  3V 

2  =  J  V  h 


J 


dw 

dz 


+  (B0-Vl)V2hQ.  =  0  at  z  =  1 


(17) 


Because  the  theory  of  exchange  volatility  also  holds  true  for  the  present  configuration,  the  time  derivatives  would 
also  be  removed  from  Eqs  (10)  and  (11)  respectively.  Then  the  dependent  variables  will  be  expanded  in  normal  mode  as 

(w,T,S)  =[w(z),0(z),4>(z)]  exp [i(lx+my)~\  (18) 

Then,  we  get  the  following  differential  equations  for  eigen  value  problem 

/  ( D2  -  a2  J 2 IV  +  2Df  ( D2  -  a2 )  DW  +  D2f  ( D2  +  a2 )  =  Ra2G  +  —  Rs  a2<$>  (19) 

Le 
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(D2-a2)0  +  DF(D2-a2)<S>  =  -W\\-Ns(\-2z)\  (20) 

Sr(D2-a2)0  +  Le(D2-a2)<5>  =  -W  (21) 

The  linearized  boundary  conditions  are: 

W  =  DO  +  Bi[0-(l  +  Ns)z]  =  0  at  z  =  1  (22) 

f(D2+a2)w  +  Ma2[0-(\  +  Ns)z]  =  Q  at  z  =  l  (23) 

f  Cr(D2  -3a2)  DW  =  (BQ  +  a2)a2Z  at  z  =  l  (24) 

W  =  0,  DW  =  0  and  DO  =  0  at  z=0  (25) 

4.  ANALYTICAL  SOLUTION 


Equations  (19)-(21)  can  be  solved  analytically  using  regular  perturbation  technique  (  RPT ) ,  subjectcd  to  boundary 
conditions,  equations  (22)-(25).  The  variables  are  expressed  in  terms  of  the  small  wave  number  to  study  the  validity. 


{W,0,^)  =  Yi(a2)i  (26) 

1=0 

Substitution  of  Eq.(26)  into  Eqs.(19)-(21)  and  the  boundary  conditions  (22)  -  (25) 

/  D4W0  +  2 Df  D3W0  +  D2]  D2W0  =  0  (26) 

D2Oo  +  DFD2<t>o=-f(z)W0  (27) 

Sr  D200  +  Le  D2<&0  =  -W0  (28) 

where 

f(z)  =  [l-Ns{l-2z)]  (29) 

The  required  boundary  conditions  (22)  -  (25),  takes  the  form 

W0  =  DW0  =  D®0  =0  at  z  =  0  (30) 

f(l)D2W0  =  D®0  =0  at  z  =  1  (31) 

Then  solutions  for  equations  above  are 

©0  =  1  and  W0  =  0  (32) 

First-  order  equations  (19)  -  (21)  become 
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r 


D*WX  +  2 B  D%  +  B 2  D2Wx  = 


1 


■\ 


R  h - Rs 

v  Le  J 


Exp[-B{z- 1/2)] 


D29l=\  +  DF-f{z)Wx 


fl  =  Sr  +  Le  -  Wx 


The  boundary  conditions  (22)  -  (25) 


become 


WX  =  DWX  =0  at  z  =  0 


f(\)D2W0+Ma2[Q-(\  +  Ns)z[  =  0  at  z  =  1 


~f(\)D2Wx-^Z0=Oat  z  =  1 
Cr 


The  general  solution  of  (33)  is 

f  1  ^ 


w1  = 


i?H - Rs 

Le 


Cx  +  C2z  +  C3e~Bz  +  C4  ze-Bz  +  Exp[-B  ( z  - 1/2)] 

2  B~ 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


(39) 


where 


e~B/2 

(  22e~B  +B-1  ( 

p~B/2 

r  — 

(  2-2\eB  +3B  +  \\BeB( 

2  B2 

[  1 0  -  5e~B  +  5e~2B  -  7  J 

’  2  -2  B2 

[  7  e2B-B3eB-B  ) 

eB/2 

f  e~B-B2- 1  ] 

PB/2 

r  - 

(  2-2\eB +\2 B3  +  B2  ^ 

2  B3 

[\  +  4e2B  -4eB  -B2eB  ) 

’  4  2  B2 

[22eB  -3e2B  +  5 B2eB  -\ ) 

The  differential  Equations  (34)  and  (35)  involving  D2dx  and  I)  2<T>  gives  the  solvability  condition  which  is 

given  by 


l 

J(l  +  f(z))Wx  dz  =  (1  +  Df  +  Sr+  Le)  (40) 

o 

The  expressions  for  Wx  is  back  substituted  into  Eq.  (40)  and  integrated  to  yield  an  relation  for  the  critical 
Rayleigh  number  and  Marangoni  number  are  obtained. 

5.  RESULTS  AND  DISCUSSIONS 

In  this  paper,  we  focused  on  the  combined  effects  of  the  viscosity  parameter  and  the  internal  heat  production  in  the 
presence  of  thermo  and  thermal  diffusion  parameter  on  Rayleigh-Benard  convection  (due  to  buoyancy)  and  Marangoni 
convection  (due  to  surface  tension).  A  linear  stability  study  was  performed,  the  boundaries  are  considered  rigid,  but 
permeable,  and  insulated  to  temperature  disturbances,  and  the  obtained  eigenvalue  problem  was  solved  using  the 
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perturbation  procedure.  In  carrying  out  the  analytical  computations, 

5.1  Rayleigh-Benard  Convection 

Internal  heat  generation,  Ns  and  variable  viscosity  viscosity,  B  effects  is  investigated  theoretically  using  perturbation 
procedure,  critical  Rayleigh  number  is  obtained  and  the  marginal  stability  curves  in  the  plane  is  identified.  Marginal 
stability  curves  will  shift  upwards  when  the  value  of  Ns  increases  and  curves  will  shift  downwards  when  the  value  of 
B  increases  which  state  that  the  internal  heating  in  the  fluid  layer  stabilizes  the  system  and  the  temperature-dependent 
viscosity  destabilizes  the  system.  From  figures  2  and  3,  the  trends  of  stability  for  two  parameters  which  are  the  thermal 

diffusion  parameter,  Sr  and  thermo  diffusion  parameter,  Df  that  exist  in  a  double  diffusive  is  investigated.  The  stability 

curves  for  the  effects  are  shown  in  figure  2  and  figure  3.  As  can  be  seen  clearly  in  figure  2,  R  increases  when  Sr 
increases  and  thus  enhance  the  onset  of  convection  due  to  the  decrease  of  temperature  flux.  Meanwhile,  for  Dufour 
parameter,  it  shows  that  an  increase  of  Dufour  parameter,  Df  will  increase  the  R(  as  illustrated  in  figure  3. 


Figure  2:  Effects  of  B,  Ns  &  DF  to  Critical  Rayleigh  Number  with  Sr  =  0.5. 


Figure  3:  Effects  of  B,  Ns  &  DF  to  Critical  Rayleigh  Number  with  Sr  =  1 . 
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Figure  4:  Effects  of  B,  Ns  &  DF  to  Critical  Marangoni  Number  with  Sr  =  0.5. 


Figure  5:  Effects  of  B,  Ns  &  D F  to  Critical  Marangoni  Number  with  Sr  =  I . 
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Figure  6:  Effects  of  B ,  Ns  &  Cr  to  Critical  Marangoni  Number  with  Sr  =  I  &  DF  =  1 . 

5.2  Marangoni  Convection 

Marangoni  convection  (due  to  surface  tension)  on  the  onset  of  steady  convection  in  a  double-diffusive  tluid  layer  is 
investigated  theoretically.  Variable  viscosity  parameter,  B  and  internal  heat  generation,  Ns  effects  with  of  Soret  and 
Dufour  effects  are  tested,  the  critical  M (.  increases  as  shown  in  figure  4  and  5.  The  trend  of  stability  for  the  Soret  effect, 

Sr ,  internal  heating  source,  Ns  and  Dufour  effect,  Df  is  also  shown  in  the  same  figures.  The  variation  of  increasing  the 
Soret  and  Dufour  effects  increases  the  critical  Marangoni  number.  This  observation  reveals  that  in  Marangoni  convection, 
variable  viscosity  parameter,  B  is  to  accelerate  the  arriving  of  convection  and  destabilize  the  system.  Meanwhile,  Ns , 
Sr  and  Df  delays  the  arriving  of  convection  and  stabilize  the  system.  The  effect  of  Crispation  number,  Cr  on  the 
marginal  stability  curves  will  shift  downwards  when  the  value  of  B  increases  and  an  increase  in  value  of  Cr  is  to 
decrease  the  value  of  M c  and  thus  making  system  more  unstable  as  illustrated  in  figure  6. 


Figure7:  Critical  R  Versus  Lefor  Different  Values  of 

B  with  Sr  =  1  =  DF . 
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Figure  8:  Critical  Rc  Yersus  Rs  for  Different  Values  of 

B  with  Sr  =  1  =  DF . 

5.3  Lewis  and  Solutal  Rayleigh  Effect. 

Solutal  Rayleigh,  Rs  effect  and  Lewis  number,  Le  due  to  the  solutal  balances  in  a  binary  tluid  where  Le  is  the  ratio 
between  characteristic  diffusion  of  mass  and  lengths  for  diffusion  of  heat  and  Rs  is  the  Rayleigh  number  representing  the 
solutal  balances.  The  effects  to  the  system  were  shown  in  figure  7  and  figure  8.  Figure  7  represents  the  Le  effect  on  the 
critical  Rayleigh  number,  Rc  against  Ns  and  B  .  It  is  observed  that  Le  shows  the  same  results  as  Ns  where  an  increase 
of  the  effects  values  make  the  system  stabilized.  The  increasing  values  of  Rs  ,Ns  and  B  is  shown  together  in  figure  8. 
An  increase  of  Rs  make  the  system  more  stable.  From  all  the  figures  shown,  the  system  is  more  stable  when  Rs  ,Ns  and 
Le  higher  values  and  the  system  is  more  unstable  for  higher  values  of  B. 

6.  CONCLUSIONS 

The  combined  effect  of  internal  heat  generation  and  variable  viscosity  parameter  in  the  presence  of  thermo  and  thermal 
diffusion  parameters  is  being  examined  in  this  research  paper  where  results  show  that  an  increasing  value  of  B  accelerates 
the  onset  of  Rayleigh-Benard  convection  and  opposite  effect  are  shown  case  of  Marangoni  for  small  values  of  B .  The 
effects  Rs  ,Ns,LeanA  Df  will  decelerates  the  onset  of  convection.  Rs  ,  Ns,  Le and  Dr  stahi I ize  the  system  in  both 
type  of  convection(  Rayleigh-Benard  and  Marangoni). 
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